Spatiotemporal pulses in a liquid crystal optical oscillator 
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A nonlinear optical medium results by the collective orientation of liquid crystal molecules tightly 
coupled to a transparent photoconductive layer. We show that such a medium can give a large 
gain, thus, if inserted in a ring cavity, it results in an unidirectional optical oscillator. Dynamical 
regimes with many interacting modes are made possible by the wide transverse size and the high 
nonlinearity of the liquid crystals. We show the generation of spatiotemporal pulses, coming from 
the random superposition of many coexisting modes with different frequencies. 

PACS numbers: Pacs: 05.45.-a, 42.70.Df, 42.65.Sf 



Optical oscillators have been extensively studied in the 
past, in particular ring cavities with photorefractive gain 
have attracted much attention, both experimentally 
and theoretically 0, 01 . Here we present a different op- 
tical medium that results from the collective orientation 
of liquid crystal molecules tightly coupled to a photocon- 
ductive layer and pumped in a two-wave mixing config- 
uration [a, |6|. We show for the first time that, when in- 
serted in a ring cavity, the liquid crystal light valve has a 
gain large enough to overcome the losses, thus resulting 
in an unidirectional optical oscillator. The wide trans- 
verse size and the high nonlinearity of the liquid crystal 
light valve allows us to explore dynamical regimes re- 
mained inaccessible up to date, where a huge number of 
modes are interacting. By changing the control parame- 
ters, the system displays either regular regimes, such as 
the alternation of Gauss-Laguerre cavity modes and the 
emission of out of axis large order modes, or space-time 
chaotic behaviors like the generation of high amplitude 
spatiotemporal pulses appearing in random space points. 

In this letter, we focus our attention mostly on the 
spatiotemporal pulses. We show that they are confined 
along the three space directions and that they behave 
as coherent structures, with a characteristic size and 
life-time. We present a theoretical model for the liq- 
uid crystal oscillator and, for a set of parameters con- 
sistent with the experiment, we confirm the appearance 
of spatiotemporal pulses by numerical simulations. It is 
worth, here, to distinguish these spatiotemporal pulses 
from the solitary or localized structures that arise when 
a system presents bistability between a spatially homoge- 
nous and a pattern state [tI . While spatiotemporal pulses 
originate from a coherent superposition randomly tak- 
ing place within a large number of coexisting modes, 
localized structures are isolated cells of a correspond- 
ing spatial pattern. Moreover, spatiotemporal pulses are 
confined both in space and time, while localized struc- 
tures are confined only in space and can persist for in- 
definite time if not erased by a local, and large enough. 



perturbation. Spatiotemporal pulses are somehow rem- 
iniscent of the pulsating solitons predicted in the com- 
plex Ginzburg-Landau equation (CGLE) for mode-locked 
lasers Q or of the random-phase solitons found in the 
nonlinear Schroedinger equation (NLSE) for nonlinear 
lattices 0. In our case, due to the large time scale sep- 
aration between the fast decay of the field in the cavity 
(tc ~ 10"'' s) and the slow response time of the Hquid 
crystals (r ~ 10^'^ s), the model cannot be reduced to 
a single equation for the electric field, the dynamics be- 
ing slaved by the slow evolution of the refractive index. 
Nevertheless, the model contains the main ingredients 
for the creation of spatiotemporal pulses, that is, a large 
number of modes with partially uncorrelated phases and 
interacting through the nonlinear medium. Note that 
similar mechanisms are also at the basis of the collapsing 
filaments predicted in optical turbulence [lOl] . 
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FIG. 1: Experimental setup. 

The experimental setup is schematically represented in 
FiglB The cavity, that has a total cavity length L = 240 
cm, consists of three high-reflectivity dielectric mirrors 
and a lens of / = 70 cm focal length, which enhances the 
mode stability providing a nearly spherical configuration. 
This configuration also ensures the presence of different 
longitudinal modes that oscillate at different frequencies 
and thus can lead to the formation of spatiotemporal 
pulses The gain medium is a liquid crystal light 

valve with one of the walls made of a thin slice, 1 mm 
thickness, of the photoconductive Bi2SiO20 (BSO) crys- 
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tal [01 . An external voltage Vq is applied to the light 
valve by means of transparent electrodes deposited over 
the glass window and the external side of the BSO crys- 
tal. The BSO acts as a transparent photoconductor, thus 
modulating the voltage across the liquid crystals as a 
function of the intensity of the light passing through the 
cell [ij. The working point of the liquid crystals is fixed 
at Vq ~ 20 V, frequency 25 Hz. The thickness of the 
liquid crystal layer is d = 14 fim and the lateral size of 
the cell is 20x30 mm^. The cell is pumped by an en- 
larged and coUimated (10 mm diameter) beam from an 
Ar~^ laser (A — 514 nm), intensity /p ~ 2 mW/cm'^. 
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FIG. 2: a) Experimental phase diagram and instantaneous 
snapshots taken at b) c) Vb = 10 V ^ Ip — 2.4 mW/cm?, d) 
e) Vo = 14 V, Ip = 2.0 mW/cm^, f) Vo = 25 V, Ip = 1.3 
mW/cm?, g)Vo = 20 V, Ip = 2.0 mW/cm^; h) space-time 
plot obtained for the same parameter values as in g). 



The light amplification in the cavity is based on two- 
wave mixing (2WM) interactions in the liquid crystals 
0. The pump beam and the cavity axis are an angle 
9 = 2) mrad. The cavity field E,. is spontaneously gen- 
erated and has a frequency almost equal to that of the 
pump field Ep, so that the 2WM is nearly degenerate. In 
fact, we observe a frequency detuning of a few Hz, which 
is automatically selected by the cavity in order to maxi- 
mize the gain and that corresponds to a refractive index 
grating slowly moving inside the liquid crystals. It is well 
known that a moving grating enhances the two- wave mix- 
ing gain for photorefractive crystals [l^ , however this fre- 
quency shift was not observed before for a liquid crystal 
light valve. Also note that the liquid crystal light valve is 
a thin medium, so that the 2WM here takes place in the 
Raman- Nath regime [13] ■ The Fresnel number of the cav- 
ity, which is the ratio of the area of the diffraction limiting 
aperture to the area of the fundamental Gaussian mode, 



is controlled by a diaphragm placed in front of the light 
valve and can be changed from _F = 1 to approximately 
F — 500, which implies changing from a single trans- 
verse mode oscillation to a regime where a huge number 
of modes are interacting. An important feature, that 
differentiates our system from previous photorefractive 
cavities, is that changing the voltage Vq changes the uni- 
form refractive index ric of the liquid crystals, and thus 
the frequency detuning between the lowest order cavity 
mode and the pump beam. In order to compensate the 
detuning, the cavity adjusts its emission by changing the 
length of the out of axis transverse wave vectors [14| , so 
that the number of active modes can be changed by keep- 
ing the Fresnel number fixed and by varying the voltage 

For the purpose of visualization, a small fraction (4%) 
of the cavity field is extracted by a beam sampler and di- 
rected to a CCD camera. A lens is inserted on the path 
of the extracted beam in such a way that the transverse 
intensity distributions recorded by the CCD are equiva- 
lent to the intensity distributions at the entrance plane of 
the liquid crystal light valve, z = 0, by taking the z coor- 
dinate along the cavity axis. For a fixed Fresnel number 
F — 500, and by changing Vq and the pump intensity /p, 
we have determined the experimental phase diagram, as 
reported in Figl2^. Cavity mode oscillations are in the 
larger grey area whereas spatiotemporal pulses are in the 
darker area. The low Vq regimes are similar to those pre- 
viously reported for low F photorefractive cavities, with 
the alternation of low order Gauss-Laguerre modes [H, 0] . 
The transition to the high Vq regimes is accompanied by 
the emission of high order and out of axis symmetrical 
modes. For intermediate values of Vq and Ip, a large 
number of modes interact through the nonlinear medium 
and give rise to the formation of spatiotemporal pulses, 
appearing as large intensity peaks over a lower amplitude 
and "speckle-like" background. Instantaneous snapshots 
of the transverse intensity distributions Ic{x,y), where 
X, y are the coordinates in the transverse plane, are dis- 
played in Figl2jD-g. The transverse size of the oscillating 
field increases with Vq, up to a large ring for high Vq. 
The spatiotemporal pulses appear when a large number 
of modes is populating the whole size of the area illu- 
minated by the pump beam. A typical snapshot cor- 
responding to this case is displayed in Fig[2^ and the 
corresponding spatio-temporal plot is shown in FiglJh. 

To investigate the dynamics of spatiotemporal pulses 
we have fixed Vq = 20.3 V r.m.s. and Ip = 2.0 mW/cm^ 
and we have recorded several movies. The spatiotem- 
poral pulses are identified by applying on each frame a 
threshold of for instance 3 times the average intensity, 
calculated as Ic = Y^k=i^c{xk,yk)/N, with Ic{xk,yk) 
the intensity of the cavity field in the fc-th pixel point 
and N the total number of pixels composing the image. 
The extension of the pulses in the z direction is investi- 
gated in the following way: we have divided in three parts 
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the beam extracted from the cavity, directed each part 
through an optical path of different length and recorded 
the three transverse intensity distributions by using sepa- 
rate CCDs. The electronic delay time of the trigger being 
negligible with respect to the liquid crystal response time, 
we can consider that the CCDs simultaneously record the 
slaved field evolution in three different z positions. By 
inspecting the movies, we have found that the typical 
longitudinal extension of the pulses is approximately 30 
cm, so that we have chosen to record the intensity distri- 
butions Ic{x,y,Zi) i = 1,2,3 at distances corresponding 
to zi = 0, Z2 — 5 and 2:3 = 32 cm. We cut the spatio- 
temporal plots along the x,y ov t direction and we obtain 
the ID spatial or temporal profile of the pulses. The spa- 
tial, temporal, profiles of a pulse are displayed in FiglSHH 
respectively. It can be seen from FigOJ; and Fig|3|D that 
at the distance z^ the pulse has disappeared. By taking 
the half height width of the pulses with Jc(x, y, z) > 3/c 
and by averaging over more than one hundred profiles, we 
find that the typical transverse size of a pulse is 250 ± 50 
fim whereas its average lifetime is around 0.5 ± 0.1 s. 
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FIG. 3: Spatial profiles recorded in a) 21, b) 22 and c) 2:3. 
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FIG. 4; Temporal profiles recorded in a) 22 and b) 23. 

The model, derived by coupling the Maxwell equa- 
tions for the cavity field with a Debye relaxation equa- 
tion for the refractive index [l5| . takes into account the 
Kerr nonlinearity of the medium as well as the two- 



wave-mixing mechanism of photon injection inside the 
cavity. The liquid crystal light valve is positioned in 
z — 0, perpendicularly to the cavity axis and f± de- 
notes the coordinates in the transverse {x,y) plane. The 
refractive index n(r±,t) satisfies the equation rdtn — 
nc — n + lQ'\/']_n — a\E{f±,t)\'^, where a > is the nonlin- 
ear coefficient of the light valve, ric is a constant value de- 
termined by the voltage Vq, E{rj_,t) = _E;pe*('^^'-^-"p*' + 
^^g4(fe^z-wpt) _|_ ^j^g |;o^al electric field, with Ep the 

pump amplitude and E^, = Ec{r±, t) the complex ampli- 
tude of the cavity field. The intensity, \E\^ = lEp]"^ + 

|i?cP + {E*Ec/2)e~'-^^'^^ + c.c, gives rise to a refractive 
index change with two components, one varying slowly in 
space, the other corresponding to a spatial grating with 
wave number k \ — k„ — kf We thus write the refractive 



index as n = 



- a\Ep\^ + no + (7ii/2)e 



+ C.C, 



By substituting this expression in the equation for n we 
obtain 



= (-l + ^o'Vi)n„-a|£;,|^ 

~dt 



(1) 



= [-l + lt^i+Qki_Y)ni-2il'^k^-^ni 



- aE;E, 

for the slowly varying fields np e ni. The wave equation 
for the cavity field is written by adopting, as usual, the 
slowly varying amplitude approximation and by consid- 
ering a planar cavity. Due to the large scale separation 
between the medium response time and the cavity round- 
trip time, we neglect the time derivative and obtain 



BE, 
dz 



2k c 

ikcniW{z)Ep 



ikcHoW (z) 



Ec 



(2) 



where W{z) = 1 in the liquid crystal light valve, z = to 
d, and elsewhere, d is the frequency detuning between 
the lowest order cavity mode and the pump field and 7c 
is the cavity loss rate. When deriving this equation we 
have considered rig and rii small. 

Inside the liquid crystal light valve we neglect the 
transverse Laplacian and we integrate in z the field equa- 
tion, Eq.([2]), thus obtaining 

Ec{d) = e*'==*'« [Jo{2kcd\ni\)Ec{0) + iJi{2kcd\ni\)Ep{0)] , 

(3) 

which accounts for Raman-Nath diffraction [l3| of the 
pump in the 2WM process. Outside the liquid crystals 
the field evolution is governed by difi^raction, thus the 
transverse Laplacian has to be retained. By considering 
that the cavity field has to satisfy the periodic boundary 
conditions imposed by the cavity, and by taking into ac- 
count the presence of the lens, it can be shown [16] that 
the field at the entrance of the cell is given by 
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(4) 



m=l 



where 



B 



and 
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s^'^'="Vo(2M|ni|) 
(5) 



C 



-Vi +i 



^e^'^'="«Ji(2fc,d|ni|), 
(6) 

with Li the distance between the hquid crystal hght valve 
and the lens, L2 — L — Li. The interpretation of the 
above equations is quite simple. The first term of the 
sum, E^°'^ = iCEp{0), is Eq.® with Ec{0) = plus the 
field evolution in the cavity, thus accounts for the field 
generated by the pump in one cavity round-trip. The 
following terms sum up the field evolutions on the suc- 
cessive cavity loops. Inside the liquid crystal light valve, 
the operator B gives the same field equation as Eq.Q 
but with Ep(0) = 0. At variance with usual treatments, 
where the mean-field approximation is used to eliminate 
the z dependence of the field, our model keeps this depen- 
dence, which is necessary to accounts for the formation of 
3D patterns. Previously, the z behavior was considered 



for optical parametric oscillators [17[ and in dispersive 
self- focusing media ^131, where optical bullets have been 
predicted as the 3D analogous of optical solitons. 




FIG. 5: Numerical pulse profile along z. 
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FIG. 6: Numerically calculated field distribution in the x — z 
plane at a fixed time; the dashed line marks the lens position 
(Li = 1.3 m, L2 = 1.1 m and / = 0.7 m). 

We have performed numerical simulations of the model 
equations, Eq.([T]) and Eq.(|3|) in the 2D+1 approximation. 



that is, by keeping the time and the x,z spatial depen- 
dences. The parameters are chosen from the typical val- 
ues of the experiment, t — AO ms, Iq = 30 fj-m,lp — 2 
mW/cm^ and a = 4 crri^/W. The number of terms in 
the sum is truncated to the number of round-trips given 
by the average lifetime of photons in the cavity. In Fig[5] 
is shown the profile of a spatiotemporal pulse along z. 
In Figl6]we show the intensity distribution calculated in 
the X, z plane at a fixed time. In agreement with the 
experiment, spatiotemporal pulses appears as filaments 
confined both in the transverse direction and along the 
light propagation direction z. 

In conclusion, we have shown a new type of nonlin- 
ear optical oscillator, which includes a thin liquid crystal 
light valve as the gain medium, and we have given evi- 
dence of the formation of spatiotemporal pulses. 
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